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THE RELATIONSHIP BETWEEN A STRIP WIENER HOPF 
PROBLEM AND A LINE RIEMANN HILBERT PROBLEM 

ANASTASIA V. KISIL 


Abstract. In this paper the Wiener-Hopf factorisation problem is presented 
in a unified framework with the Riemann-Hilbert factorisation. This allows 
to establish the exact relationship between the two types of factorisation. In 
particular, in the Wiener-Hopf problem one assumes more regularity than for 
the Riemann-Hilbert problem. It is shown that Wiener-Hopf factorisation can 
be obtained using Riemann-Hilbert factorisation on certain lines. 


1. Introduction 


The Wiener-Hopf and the Riemann- Hilbert problems are a subject of many 
books and articles [2S1IS112] • The similarities of the two techniques are easily visible 
and have been noted in many places. Nevertheless, to the author’s knowledge there 
was no systematic study of the exact relationship of the two methods. To fill this 
gap is the purpose of this article. 

It has been suggested in m Chapter 4.2] that the Wiener Hopf equation are 
a special case of a Riemann-Hilbert equation. Specifically, the Riemann Hilbert 
problem connects boundary values of two analytic functions on a contour and the 
Wiener Hopf equation is defined on the strip of common analyticity of two func¬ 
tions. In the simplest case, both methods use the key concept of functions analytic 
in half-planes. The additional regularity for the Wiener-Hopf equation allows to 
express the solution in more simple terms than the Riemann-Hilbert equation. To 
be well-defined the Riemann-Hilbert problem requires some additional regularity, 
e.g. the coefficients need to be Holder continuous on the contour. 

In a different book [B] Chapter 14.4] it has been stated that the Wiener-Hopf 
equations results from a bad choice of functions spaces and instead a Riemann- 
Hilbert equations should be considered. Confusingly, those Riemann-Hilbert equa¬ 
tions are sometimes referred to as a Wiener-Hopf equations. Historically, there has 
been insufficient interaction between the communities using the Wiener-Hopf and 
the Riemann-Hilbert methods, this results in obscuring disagreements in termino¬ 
logy and notations. Such differences if not reconciled properly have a tendency to 
widen. 

We consider the following problem as an illustration. Given a function F(t) on 
the real axis: 


m 


It 2 - {\i + l)t - \i + | (t-i + 2) 


t 2 — + 1 


it + i + f) ’ 


( 1 . 1 ) 


find two factors F+(t) and F (t), which have analytic extensions in the upper and 
lower half-planes respectively. In this rare case the factorisation can be obtained 
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by inspection 
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F+(t)F-(t). 
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(t — i + 2) 


( 1 . 2 ) 


These functions are depicted in Figure [Q We will comment on this example in both 
(the Riemann Hilbert and Wiener- Hopf) frameworks at the end of this paper. 

The structure of the paper is as follows. Section [2] collects some preliminaries. 
In Section [3l the usual theory of the Wiener-Hopf equation is recalled. It is presen¬ 
ted to highlight the differences with the theory of the Riemann-Hilbert equations, 
Section 0] Section [5] describes at the relationship between the two methods in 
the context of integral equations. In Section E] the solution to Wiener-Hopf and 
Riemann-Hilbert factorisation is re-expressed in terms of the Fourier transforms 
instead of the Cauchy type integrals. This allows to prove theorems about the 
exact relationship of Wiener Hopf and Riemann- Hilbert factorisation. 


2. Preliminaries 

This section will review important properties of the Fourier transform, the 
Cauchy type integral and the relationship between the two. We also states the 
generalised version of Liouville’s Theorem together with analytic continuation, in 
the form which will be used later. 

The Fourier transform of a function f(t) is defined as 

OO 

F( x) = -L= I f(t)e lxt dt, —oo < x < oo. (2.1) 

V27T J 

— OO 

Throughout this paper the Fourier image and original function will be denoted by 
the same letter but they will be upper and lower case respectively. The inverse 
Fourier transform is given: 

OO 

f(t) = f F{x)e~ lxt dt , — oo < t < oo. 

v27t J 

— OO 


It is an important fact that if f(t) 6 /^(R) then F(t) 6 L^M.). This makes the 
space of square integrable functions very convenient to work in. Moreover, Fourier 
transform is an isometry of L 2 (M) due to Plancherel’s theorem. 

The Fourier transform need not be confined to the real axis, as long as the integ¬ 
ral (12.111 is absolutely convergent for a complex x. On an open domain consisting 
of such parameters x, the Fourier transform is an analytic function. The Paley- 
Wiener theorem (see (12.61) (12.81) below) gives the conditions on the decay at infinity 
of f(t) to ensure such analyticity in different strips and half-planes. 

Another key concept is the Cauchy type integral. Let F(r) be integrable on a 
simple Jordan curve F, the integral 


1 

2ni 


F(t) 

T — Z 


dr, 


is called a Cauchy type integral and defines an analytic function on the complement 
of L. If L divides the complex plane into two disjoint open components then it 
makes sense to consider two functions F~ and F + on the respective domains. In 
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Figure 1. A full colour image of functions (11.21) . The top picture 
shows F(z) followed by F + (z) and F~(z). We use a colour scheme 
developed by John Richardson. Red is real, blue is positive ima¬ 
ginary, green is negative imaginary, black is small magnitude and 
white is large magnitude. Branch cuts appear as colour discon¬ 
tinuities. Produced using MATLAB package zviz.m. 
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particular, for the real line we use the notation: 

OO 


1 

2iri 



F+( Z ) 

F~(z) 


if Imz > 0, 
if Imz < 0. 


( 2 . 2 ) 


— OO 

The relationship of the Cauchy type integral to the Fourier transform is outlined 
below. 


OO 



— OO 
OO 



— OO 


dr = -L [ f(t)e izt dt , 

t-z y/2jr J W 
0 


0 



— oo 


f{t)e izt dt , 


if Imz >0, (2.3) 

if Imz < 0. (2.4) 


The above integrals are the two ways of arriving at functions analytic in half-planes. 

It was already mentioned that L 2 (K) is a very convenient with regards to the 
Fourier transform. The Holder continuous functions produce well-defined boundary 
values of the Cauchy integral. Thus, their intersection [6] § 1.2]: 

{{0}} = L 2 (K) H Holder, 


turns out to be very useful for both the Wiener-Hopf and the Riemann-Hilbert 
problems. The pre-image of {{0}} under the Fourier transform is denoted {0}. 

Formulae (12.31) (12.41) show the significance of functions which are zero on the 
positive or the negative half lines. Given a function on the real axis we can define 
the splitting 



if t > 0, 
if t < 0, 


/-(*) 


0 if t > 0, 
-/(f) if t < 0. 


(2.5) 


We will say that if / G {0} then f+(t) G {0, oo} and /-(f) G {—oo, 0}. 

In the rest of the paper we will need to refer to functions that are analytic on 
strips or (shifted) half-planes. Following [6] § 13], we define / G {a} if e~ ax f € {0}, 
that is a shift in the Fourier space. Finally, / = /+ + /_ G {a, 6} if /+ G {a} and 
/_ G {6}. From the definition of /+ and /_ it is clear that also f + G {a, oo} and 
/_ G {-oo,6}. 

The Fourier transform of functions in the class {a, 6} is denoted {{a, 6}}. The 
celebrated Paley-Wiener theorem states that the following is equivalent : 


F + (z) analytic in Imz > a 
F-(z) analytic in Imz < b 
F(z) analytic in a < Imz < 6 
Recall, a convolution of two functions 


OO 



<(=> F+(z) G {{a,oo}}, 

(2.6) 

F-{z) e {{-°o> & }}, 

(2.7) 

<=► F(z)e{{a,b}}- 

(2.8) 

the real line is given by 


t — s)g(s) ds. 

(2.9) 


—OO 


One of the key properties, which make Fourier techniques useful in integral equa¬ 
tions, is that the Fourier transform of the convolution is the product of the Fourier 
transform of the functions /(t) and g(t), i.e. H{x) = F(x)G(x). 

If g G {a, 6} and / G {a,/3} (a < b, a < /3) then h(x) will be in the class 
{max(a, a), min(6, /?)} [61 § 1.3] provided 

max(a, a) < min(6, /3). 

After taking the Fourier transform this will become: 

H(z) = F(z)G(z), 
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and H(z ) £ {{max(a, a), min(6, /3)}} i.e analytic in the strip 

max(a, a) < Im {z) < min(6, /3). (2-10) 

We will need a generalised version of Liouville’s Theorem together with analytic 
continuation: 


Theorem 2.1 ([B] § 3.1]). If functions F\{z), ^ 2 ( 2 ) arc analytic in the upper and 
lower half-planes respectively with exception of zq = 00 , Zk (k = 1,2... n) where 
they have poles with principal parts: 


G 0 (z) 



c°z- 


-0 > 



c 


k 

m k 


{z - z k ) m >° ’ 


( 2 . 11 ) 

( 2 . 12 ) 


with Fi(z) and F^z) equal on the real axis, then they define a rational function on 
the whole plane: 

F(z) = c + Gq(z ) + Gfe ^ , 

where c is an arbitrary constant. Poles Zk can lie anywhere on the half-planes or 
on the real axis. 


3 . The Wiener-Hopf Method 

In this section the Wiener- Hopf method is presented in the way it appears in 
most papers. This will be revisited to establish the relationship with the Riemann- 
Hilbert method. For some applications of Wiener-Hopf in areas like elasticity, crack 
propagation and acoustics see mmmm- 

The next theorem provides a constructive existence theorem for the additive 
Wiener Hopf decomposition (in terms of a Cauchy type integral). 

Theorem 3.1 f [HI Ch. 1.3]). Let f(a) be a function of variable a = a + it, 
analytic in the strip r_ < r < t + , such that ffcr + ir) < C'|cr| _p , p > 0 as \a\ —> 00 , 
the inequality holding uniformly for all r in the strip r_ + e < r < r + — e, e>0. 
Then for t_ <c<r<d<r + , 


/(a) = /-(a) +/+(a), 


00 -\-id oo+zc 

/-(«) = ; /+<«) = f- ! P^-di, ( 3 . 1 ) 


27ti J f — a 

—00 -\-id 

where /_ is analytic in t < t + and /+ is analytic in r > . 


27ri J f — a 

— 00+2C 


We wish to highlight the simplicity of demonstration of this result. Indeed, to 
prove EH) one applies Cauchy’s integral theorem to the rectangle with vertices 
±a + ic , ±a + id. From the assumption as regards to the behaviour of f(a) as 
|<t| —> 00 in the strip, the integral on er = ±a tends to zero as a —^ 00 and the result 
follows. 

The next theorem is a useful variation of the previous theorem, obtained by 
taking logarithms. This provides a way to achieve the multiplicative Wiener-Hopf 
factorisation. 


Theorem 3.2. [H] Ch. 1.3] If log K satisfies the conditions for theorem, 1,9.11 in 

particular that K(a) is analytic and non-zero in the strip and K{a) —> 1 as \cr\ —> 00 
then K{a) = K + (a)K-(a) and K + , K_ are analytic, bounded and non-zero when 
t > T—, t < t_|_ respectively. 
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The above factors are unique up to a constant HU. In other words, if there are 
two such factorisations K = K + K_ and K = P + P_, then: 

K + = cP+ and /<'_ = c _1 P_, 

where c is some complex constant. This can be seen by applying analytic continu¬ 
ation to I\ + /P + and P_/AT _ and then using the extended Liouvillc theorem 12.11 


4. The Riemann-Hilbert Method 

In this section the main differences from the Wiener-Hopf method are presen¬ 
ted. They are rooted in the conditions imposed on the functions and the resulting 
formula. For applications of the Riemann- Hilbert problem the interested reader is 
refereed to W\- 

The Riemann-Hilbert Problem can be stated as follows: 

Problem 4.1 (Riemann- Hilbert). On the real line two functions are given, H(t) and 
non-zero D(t) with 

D(t) - 1 e {{0}} and H(t) £ {{0}}. 

It is required to find two functions F± (z) analytic in the upper and lower half-planes 
such that their boundary values F ± (t) on the real line satisfy two conditions: 

(1) P ± (t) belong to the classes {{0,oo}} and {{—oo, 0}}; 

(2) The identity holds: 

D(t)F~(t) + H(t) = F + {t), for all real t. (4.1) 

As it has been seen for the Wiener-Hopf method the key steps in the solution is 
the additive splitting or jump problem (Thm. 13.11) and the factorisation (Thm. EOll . 
First splitting is addressed as before: 

Problem 4.2 (Jump problem). On the real line a function F(t) £ {{0}} is given. 
It is required to find two functions F ± (z) analytic in the upper and lower half¬ 
planes with boundary functions on the real line belonging to the classes {{0, oo}} 
and {{—oo,0}} and satisfying: 

F(t) = F+(t) + F~(t), 

on the real line. 


The solution is offered by the Sokhotskyi-Plemclj formula: 

OO 

x) — F~(x) = F(x), F + (x) + F~(x) = — / —— dr. 

7 n J r — x 

—oo 



(4.2) 


We shall note, that the derivation of m is more involved than the proof of 
Thm. 13.11 

Next the factorisation problem is examined. The index of a continuous non-zero 
function K(t) on the real line is: 

ind(AT(f)) = ^-( lim argA'(f)— lim arg K(t)). (4.3) 

In other words the index is the winding number of the curve (Re K(t),ImK(t)) 
l€l. Note that ind = 1. Given a function K(t ) with index k one can reduce 
it to zero index by considering a new function 

Mt) = Kit) (£) ” ■ 

The reason it is important to consider functions K{t ) with zero index is to ensure 
In K(t) is single-valued. For the rest of this paper we will assume that all functions 
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have zero index. Taking logarithms and applying the Sokhotskyi-Plemelj formula 
we get a solution to the factorisation problem |B]: 

Problem 4.3 (Riemann-Hilbert Factorisation). Let a non-zero function K(t ), such 
that K(t) — 1 £ {{0}} and ind AT(t) = 0, be given. It is required to find two functions 
K±(z) analytic in the upper and lower half-planes with boundary functions K±(t) 
on the real line belonging to the classes {{0,oo}} and {{—oo,0}} and satisfying: 

K{t) = K + {t)K~{t), 

on the real line. 

5. Relationship Between Wiener-Hopf and Riemann-Hilbert via 

Integral Equations 

This section demonstrates how one type of integral equation can either lead 
to the Wiener-Hopf problem or the Riemann-Hilbert problem, depending on the 
class of function where the solution is sought. Integral equations have historically 
motivated the introduction of the Wiener-Hopf equation [12] . 

In applications, a time-invariant process can be modelled by an integral equations 
with convolution on the half-line: 

OO 

J k(x - y)f(y) dy = g(x ), 0 < x < oo. (5.1) 

o 

Here the kernel k(x — y) £ (a, b} represents the process, g{x) is a given output and 
f(y ) is an input to be determined. To solve the equation (15.11) we complement the 
domain of x: 

OO 

J k(x - y)f(y)dy = h(x ), -oo < x < 0, (5.2) 

o 

where h(x) is unknown. Then, by applying the Fourier transform we get the equa¬ 
tion: 

F + (a)K(a)-G+(a) = H-(a). (5.3) 

Now it is time to examine the above equation more carefully and in particular 
clarify the analyticity regions. There will be two different cases considered: the first 
one will be a typical example from applications and the second the most general 
solution. We will see that the former will lead to a Wiener-Hopf equation and the 
latter to a Riemann-Hilbert equation. 

(1) In equation (15.31) K(a) and G+(a) are known, thus their region of regularity 
can be determined. From the maximal growth rate of f(x) as x —> +oo and 
h(x) as x —> —oo the analyticity half-planes are determined as in (12.61) 
(12.81) . For example, in the integral equation for Sommerfeld’s half-plane 
problem m Ch. 2.5], the known functions are in the following classes: 

G + (a ) £ {{acos0,oo}}, K(a) £ {{—a,a}}, 

and the unknown are in: 

F + (a) £ {{acos0,oo}}, H_(a) £ {{—oo,a}}. 

Here a and 9 are some constant. From (12.101) . K{a)F + {a) £ {{a cos 9, a}} 
and Equation (15.31) holds in the strip a cos 9 < Imcc < a . We obtained a 
Wiener-Hopf equation. 
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(2) We considered K(a) £ {{—a, a}} and the other functions will be assumed 
to belong to the largest possible class for (E3D to be solvable. If F + (a) £ 
{{6, oo}} the convolution will exist for b < a, so we will take equality as the 
minimal condition on regularity. Then, the maximal class [6] , in which the 
integral equation has a solution, is: 

G +(a) S {{a,oo}}, F + (a) £ {{a, oo}}, H_(a) £ {{-oo, a}}. 

Furthermore, from (12.101) we have K(a)F + (a) £ {{a,a}}. Hence, (15.dll is 
only valid on a line Im a = a and it is a Riemann- Hilbert equation. 

The above shows that the Wiener Hopf equation is the result of a better regu¬ 
larity of the function at infinity than is minimally needed for a solution to exist. 


6. Relationship between the Wiener-Hopf and the Riemann-Hilbert 

Equations 

To examine the relationship between the Wiener Hopf and the Riemann-Hilbert 
equations we will restate problems for the same class of functions and re-express 
the solution in terms of the Fourier integrals instead of Cauchy type integrals. 

We begin with the jump problem in the Riemann-Hilbert case: 

Theorem 6.1. The solution to Jump Problem \f.2\ can be expressed as: 

oo 0 


F + {z) = 


1 


V2tt . 


f(t)e izt dt, F~ (z) = — 


V2tt 


f(t)e izt dt, 


( 6 . 1 ) 


where f(t) is the inverse Fourier transform of F(t). 

Proof. The solution to the problem is given by the Sokhotskyi-Plemclj formula (HUD- 
We re-express them in terms of the Fourier integrals using 


F(x) = 


V2^ 


f(t)e ixt dt, 


and 


OO 

L S 

m J 


F(t) 


hl m 


e lxt dt + 


\[Tk 


f(t)e ixt dt , 


= F + (x) — F ( x ) 


T — X 


dr = F + (x) + F (x), 


1 ( f(t)e ixt dt - 1 






f{ty 




sgn (t)f{t)e ixt dt. 


In other words, the Sokhotskyi-Plemelj formula for the real line in terms of the 
Fourier transform is: 


F + (x ) — F (x) = F{x ), F + (x) + F (x) = 


1 




sgn (t)f{t)e ixt dt, (6.2) 


From this the result follows. 


□ 
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Remark 6.2. The problem does not change if it is shifted by real a to the classes 
F(t) 6 {{a}} with {{a, oo}} and {{—oo,a}}. 

To clarify the relation between two problems we re-state the Wiener-Hopf jump 
problem from Thm. [3711 as follows: 


Problem 6.3 (Wiener-Hopf jump problem). On a strip S = {z : a < Inr (z) < 6} a 
function F(z) € {{a, 6}} is given. It is required to find two functions F ± (z) analytic 
in the half-planes {z : a < Im (z)} and {z : Im (z) < b} respectively, which belong 
to the classes {{a, oo}} and {{—oo,6}} and satisfying: 

F(z) = F + {z) + F~(z), 


on the strip S. 


Similarly to Thm. 16.11 we find: 


Theorem 6.4. The solution of Prob. 1 6 . ,91 can be expressed as: 

oo b 

F + (z) = -j=j f(t)e izt dt F~(z) = -L j f(t)e izt dt, 

a —oo 

where f(t) is the inverse Fourier transform of F(t). 


Proof. To prove the theorem we use the solution (ED of the Riemann- Hilbert 
problem first for the line Im (z) = a and then for Im (z) = b. On the line Im (z) = a 
we obtain: 

F(x + ia) = F + (x + ia) + F~(x + ia), (6.3) 

with the factors given by: 


F + (z) = —^= f f(t)e izt dt and F~{z) 

v27T J 



f(t)e izt dt. 


(6.4) 


Note that: 


F(x + ia) — F (x + ia) = F + (x + ia), 


(6.5) 


where the left-hand side has continuous analytic extension into the strip S, thus 
the same will be true for F + (x + ia). This allows to move to the other side of the 
strip: 

F( x + ib) = F + [x + ib) + F~ (x + ib). (6.6) 

Now, an application of the Sokhotskyi-Plemelj formula gives: 


F+{z) = ^=jf{t)e izt dt. (6.7) 

b 


In other words, the analyticity of F + (z) was extended to the strip S and the result 
follows. □ 


Remark 6.5. Note, that in the above proof one could have chosen initially the 
Im (z) = b and then extended function F~(x+ib) to the line Im (z) = a. In fact, any 
line in between Im (z) = c with a < c < b could have been taken and a solution of 
the Riemann-Hilbert problem obtained. Then, the existance of analytical extension 
of both functions F + and F~ can be shown in a similar manner. 

Similarly, we describe the relationship between the factorisation in the Wiener- 
Hopf and Riemann-Hilbert setting. 
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Theorem 6.6. The solution or Prob. \f.3\ can be expressed as: 

K + (z) = exp -L= J n(t)e lzt dtj and K~{z) = exp ^ J n{t)e lzt dt 

where n(t) is the inverse Fourier transform of In K(t). 

Proof. From the ind AT(f) = 0 it follows that In K(t) is single-valued. An application 
of Thm. 16.11 to In K(t) yelds: 

oo / o \ 

li\K(t) = -^= J n(t)e' lzt dt+ I — — J n(t)e lzt dt j . 


The result follows from taking exponents of both sides of the last identity. 


□ 


To express relations between two factorisation problems we formulate the Wiener- 
Hopf factorisatiojn in a suitable form. 

Problem 6.7 (Wiener-Hopf Factorisation). A function A'(z) is non-zero on the 
whole strip S' = {z : a < Im(z) < b }, furthermore K(z ) — 1 £ {{a, b}} and 
ind K(x + ia ) = 0. It is required to find two functions Fl (z) analytic in the 
half-planes {z : a < Im(z)} and {z : Irri (z) < b} belonging to the classes {{a,oo}} 
and {{—oo, fo}} respectively, such that: 

K(z) = I\ + (z)K~ (z), 

on the strip a < Im (z) < b. 

Similarly to Thm. 16.61 we find: 

Theorem 6.8. The solution of Prob. can be expressed as: 

K + (z) = exp J n(t)e lzt dt^ and K~(z) = ex p J n(t)e lzt dt 

where n{t) is the inverse Fourier transform ofhiKft). 

Proof. The function K (z) on the line Imz = a satisfies all the assumptions of 
Thm. 16.61 thus we obtain: 

K(x + ia) = K + (x + ia)K~(x + ia), 

with: 


/ oo \ /a 

K + (z) = exp I — J K(t)e lzt dt J and K~{z) = exp I — J n(t)e 


' dt 


Since indFf(a: + ia) = 0 and K(z) zero free on the strip, it follows that ind K(x + 
is) = 0 for all a < s < b. This is again expressed as 

K(x + ia)K~{x + ia) = K + {x + ia). 

Because the left hand side has continuous analytic extension in the strip S, the 
function K + (x+ia) has the extension as well. This gives meaning to the expression: 

K(x + ib) = K + {x + ib)K~ {x + ib). 

Now, K{x + ib) satisfies all the assumptions of Thm. 16.61 thus 


K + (z) = exp 


V2tt ■ 


n(t)e izt dt . 


This provides the required factorisation. 


□ 
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Remark EH also holds here. There are some differences in the formulation of 
Thm. 13.21 and Thm. EH The most significant differenccQ is the assumption that 
ind K(x + ia) = 0. This is because Thm. EH assumes the existence of single valued 
In K (t). 

The above derivations show that the Wiener- Hopf equations is characterised by 
extra regularity, namely a domain of analyticity. Noteworthy, there are further 
applications of this feature, for example, a strip deformation. Consider a Wiener 
Hopf equation: 

H(a)$+(a) + T_(a) + C(a ) = 0. 

Assume as before that <I> + (a) and T_(a) are analytic in the upper or lower half¬ 
planes and the strip respectively. However, assume that this time A(a) and (7(a) 
have singularities in the strip . Then, by taking a subset of the strip the Wiener- 
Hopf equations can still be solved PH- 

Finally, we are going to revisit the factorisation of F(t) (11.11) . Instead of spotting 
the factors by inspection the Riemann-Hilbert formula can be used, yielding factors 
analytic in the half-planes. By inspection of singularities in the complex plane, F(z) 
has a strip of analyticity —1 + e < Im ( z) < 1/2 — e. Hence, it is also possible to 
apply the Wiener Hopf formula to obtain the factorisation. Due to the uniqueness 
of the factors and analytic continuation it follows, that both methods shall produce 
results coinciding with (11.21) and illustrated by Figure |1] 
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